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A prototype coupled Marine Atmospheric Boundary Layer 
(HABL) and Oceanic Boundary Layer (OBL) aiDdel, from physical 
models which have been separately formulated, is developed. 
The observational efforts will be directed toward coupling 
(local) point MABL and OBL features. Emphasis in the char- 
acreriza ticn of local features will be on the evolutions of 
the adjacent well-mixed layers. The approach will be to 
compare observed evolutions in the oceanic and atmospheric 
boundary layers with predictions from bulk models wherein 
the ocean influences the armosphere through the surface tem- 
perature which, in turn, is influenced by atmospheric forc- 
ing. The atmospheric forcing parameters are the local 
surface turbulent fluxes and radiation. Surface turbulent 
fluxes are controlled by the surface temperature and armos- 
pheric mixed layer parameters. Radiation is controlled by 
predicted cloud formation and dissipation. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

II. BACKGROUND 14 

A. GENERAL BOUNDARY LAYER FEATURES 14 

B. MARINE ATMOSPHERIC BOUNDARY LAYER (HABL) MODEL 23 

C. OCEANIC BOUNDARY LAYER (OBL) MODEL 31 

D. FORMULATED COUPLED BOUNDARY LAYER MODEL .... 37 

III. DATA AND MODEL RESULTS 43 

A. RESULTS FROM COUPLED MODEL 43 



1. CASE 1 (0500 1 9 May to 0500 20 Kay) . . 48 

2. CASE 2 (0500 2 0 May to 0500 21 May) . . 52 

3. CASE 3 (1900 0 4 OCT to 1 900 05 OCT) . . 57 

B. COMPARISON CF COUPLED AND ATMOSPHERIC MODEL . . 59 



IV. CONCLUSIONS/RECCMMENDATION S 62 

LIST OF REFERENCES 36 

INITIAL DISTRIBUTION LIST 38 



5 



LIST OF TABLES 



TABLE I. CouDled Model Predicted Atmospheric and 

Oceahxc Mixed Layer Quantities- Case 1. . . . 64 



TABLE II 



MABL Model Predicted Atmospheric Mixed Layer 
Quantities- Case 1 65 



TABLE III. Coupled Model Predicted Atmospheric and 

Oceanic Mixed Layer Quantities- Case 2. ... 71 



TABLE IV . 



HABL Model Predicted Atmospheric Mixed Layer 
Quantities- Case 2 72 



TABLE V. Coupled Model Predicted Aimospheric and 

Oceanic Mixed Layer Quantities- Case 3. . . . 78 



TABLE VI. 



MA3L Model Predicted Atmospheric Mixed Layer 
Quantities- Case 3 79 



6 



LIST OF FIGORES 



Figure 


1 . 


Coup''ed Environmental Model System 


13 


Figure 


2. 


Simplified Atmospheric and Oceanic Boundary 
Layer Temperature Profiles 


16 


Figure 


3. 


Mechanical Energy Budget for the Ocean Mixed 
Layer 


18 


Figure 


4 . 


Atmospheric and Oceanic Mixed Layer 
Observat"^ ons dur' ng CEWCOM-78 


21 


Figure 


5 . 


Simplified Flow Diagram of the Boundary Layer 
Model : . . . 


24 


Figure 


6 . 


Schematic of Input, Prescription and 
Computing Steps in MABL Pred*ot*on 


27 


Figure 


7. 


Schematic of Input, Prescripri on and 
Computing Steps in OBL Prediction 


38 


Figure 


8 . 


Schematic of Atmospheric and Oceanic Coupled 
Model 


42 


Figure 


9 . 


Anchorage Position of R/V ACANIA in the 
Vicinity of San Nicolas Island from 5/19/78 
to 5/2 0/7 8 


45 


Fi gure 


10. 


Positions of R/V ACANIA from 5/18/78 to 
5/25/78 


46 


Figure 


1 1 . 


PosiT.ions of R/V ACANIA from 10/4/76 to 
10/6/76 


47 


Figure 


12. 


Time Evolution of the Sea-Surface 
Temperature- Case 1 


66 


Figure 


13. 


Time Evolution of Ocean Mixed Laver Depth 
(MLD)- Case 1 1 


66 


Figure 


14. 


Selected Ocean Vertical Temperature Profiles- 
Case 1 


67 


Figure 


15. 


Selected Ocean Vertical Profiles of the 
Horizontal Wind Driven Currents- Case 1. . . 


67 


Figure 


16. 


Selected Change in Ocean Vertical Salinity 
Profiles- Case 1 *. . 


68 



7 



69 



Figure 17. Graphical Display of th$ Coupled Model 

Predicted Atmospnerrc Mixed Layer Quantities- 
Case 1 

Figure 18. Graphical Display of the HABL Model Predicted 
Atmospheric Mixed Layer Quantities- Case 1. . 

Figure 19. Time Evolution of the Sea-Surface 

Temperature- Case 2 73 

Figure 20. Time Evolution of Ocean Mixed Layer Depth 

(MLD)- Case 2 7 3 

Figure 21. Selected Ccean Vertical Temperature Profiles- 

Case 2. ...... ............. 74 

Figure 22. Selected Ocean Vertical Profiles of the 

Horizontal Wind Driven Currents- Case 2. . . 74 

Figure 23. Selected Change in Ocean Vertical Salinity 

Profiles- Cass 2 75 

Figure 24. Graphical Display of -^he Coupled Model 

Predicted Atmospheric Mixed Layer Quan'^rit ies- 
Case 2 76 

Figure 25. Graphical Display of the MABL Model Predicted 

Atmospheric Mixed Layer Quantities- Case 2. . 77 

Figure 26. Time Evolution of the Sea-Surface 

Temperature- Cass 3 80 

Figure 27. Time Evolution of Ocean Mixed Layer Depth 

(MLD)- Case 3 '. . . . 30 

Figure 28. Selected Ocean Vertical Temperature Profiles- 

Case 3 81 

Figure 29. Selected Ocean Vertical Profiles of the 

Hcrizontal wind Driven Currents- Case 3. . . 81 

Figure 30. Selected Change in Ocean Vertical Salinity 

Profiles- Case 3 ‘ . 82 

Figure 31. Graphical Display of the Coupled Model 

Preaicted Atmospneric Mixed Laver Quantities- 
Case 3 ; 83 

Figure 32. Graphical Display of the MABL Model Predicted 

Atmospheric Mixed Layer Quantities- Cass 3. . 84 



8 



Figure 33. 



Fog and Stratus Cloud Positions During Case 3 
Prediction Period 






?6 



■ 

• t * » 4. • 



IMJC*' 

*'■ ♦' • ♦ f 



,^5 




ft '^•al 

i,-7T • 




'V 



if ft’ 




? h*xdv4^ 

■i 




ACKN0WL2DS EMENT 



The author wishes to thank Professor K.L. Davidson, 
Department of Meteorology, and Dr. R.W. Sarwood, Department 
of Oceanography, Naval Postgraduate School, for their time, 
inrerest and guidence throughout this study. Dr. C.W. 
Fairall for his helpful discussions and review of this the- 
sis, and to CDR J.J. Jensen who provided motivation for this 
effort. Particular appreciation is extended tc Peter Guest 
for his assistance in programing and research. 

Finally, I would like to express my gratitude to my 
wife, Stephanie, and to my two daughters, Jamie and Lorry, 
for their patience, understanding and support during the two 
years of study. 



10 



I 



INTROD U CTI ON 



The idea that both atmospheric and subsurface environ- 
mental conditions can have a significant effect (enhancement 
or degradation) on modern weapon and sensor systems has now 
become a reality to many military leaders. This has led to 
a greater demand for more reliable and accurate atmospheric 
and oceanic prediction schemes. Prediction of the existence 
and evolution of electromagnetic ducts are important in pre- 
dicting over the horizon surface to surface detection. Pre- 
dictions of the inversion height is important where it is 
known that optical propagation is degraded due to turbulence 
and extinction. Ocean mixed layer predictions are important 
for determining convergence zone distances, optimum source 
level depth, and cutoff frequency for acoustic surface duct 
pro pagaticn. 

The ability to provide accurate forecasts requires a 
complete understanding of the dynamical processes occurring 
in each of the atmospheric and oceanic boundary layers, and 
as important, the processes associated with coupled local 
changes in the two layers. 

Very little is known about the interactive nature of the 
adjacent turbulent boundary layers. Th? objective of this 
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thesis is to formulate a model which couples two already 
existing and proven models for the respective adjacent 
well-mixed boundary layers. This coupling is achieved by 
matching the surface boundary conditions for stress, mass 
flux and heat flux (Figure 1) . 

This coupled model may then be used to examine the pro- 
cesses of coupled local oceanic and atmosoheric mixed layers 
at time and spa-^ial scales of 12-18 hours and 50-100 km, 
respectively. Significant changes in both mixed layers can 
occur at these scales. Explaining these changes requires 
coupled evaluations of exisring atmospheric and oceanic 
well-mixed layers at a point. 

At the present tine, we only have qualitative evidence 
•^hat changes in each layer are correlated. The rigorous 
specifioat ions or descriptions of the coupling will require 
considerable interpretive efforts and, perhaps, improved 
formulation of the separate models. 
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II. ba£^SS2M2 

A. GENERAL BOUNDARY LAYER FEATURES 

The Marine Atmospheric Boundary Layar (MABL) extends 
from the surface through the capping inversion. The Oceanic 
Boundary Layer (OBL) extends from the surface to the sea- 
sonal thermocline. The air-sea interface is bordered by 
oceanic and atmospheric turbulent mixed layers which effec- 
tively insulate the quasi-geos trophic oceanic and atmos- 
pheric flows below the thermocline and above the inversion, 
respectively. The primary source of the turbulence is the 
velocity (current) and buoyancy (density) gradients created 
by interaction a* the interface. The large vertical mixing 
yields nearly homogeneous (well- mixed) wind, temperature and 
humidity profiles above the surface. In the ocean, the term 
mixed layer refers to the region of vertical uniformiry in 
mean velocity and density. Although density is a function 
of both remperature and salinity, it is usually dependenr 
upon temperature except when large near-surface salinity 
changes occur due to heavy precipitation or ice melt in high 
latitudes. At the bottom of the ocean mixed layer and at 
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the top of the atmosphere mixed layer are thin transition 
regions (thermocline cr inversion) . The mixed layers inter- 
act with the upper and lower layers at these regions by 
means of turbulent forced entrainment. Typical model pro- 
files depicting the boundary layer regions for both the 
atmosphere and ocean are shown in Figure 2. 

The well-mixed nature of both the turbulent MABL and OBL 
implies features for the vertical distribution of mean val- 
ues and vertical fluxes of wind, temperature and humidity in 
the atmosphere and density and momentum in the ocean. One 
feature is that the properties which are conserved during 
mixing can be treated as being constant with height (depth) 
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solely on the surface and inversion (entrainment) fluxes if 
advection is negligible. 

Another important source of energy, in addition to tur- 
bulent kinetic energy, in both mixed layers is radiation. 
Unlike the HA3L, most of the solar radiation does not pene- 
trate the OBL. Therefore, downward turbulent heat flux in 
the OBL is as important as the upward flux. Since the long 
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Temperature Profiles. 
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wave (infrared) radiation flux is absorbed and emitted 
within the first millimeter of the sea surface, this flux 
component is essentially anorher component of the oceanic 
surface heat flux. The tonal surface heat flux is comprised 
of sensible and latent heat and the long and short wave 
solar radiation. Both the long and short wave (solar) 
fluxes into the ocean are critically dependent upon the 
existence of clouds in the atmosphere. 

In both the MABL and OBL, buoyantly-created velocity 
fluctuations and mixing have the most pronounced effect on 
mixed layer charact e ristics. Garwood (1977) states that 
buoyancy-forced vertical mixing has a more obvious and 
direct role in the mechanical energy budget, as shown in 
Figure 3. In the atmosphere, even under near neutral condi- 
tions, buoyancy forced fluctuations and mixing can be quite 
large. These buoyancy driven energetic eddies fill the 
atmospheric boundary layer from the surface to the inversion 
where they entrain warm, dry air and bring more momentum 
into the mixed layer. If entrainment causes the HABL to 
extend above the lifting condensation level, clouds or fog 
will form within the layer. 
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Figure 3. Mechanical Energy Budget for the Ocean Mixed 
Layer . 

It is the r=-alization that buoyancy driven entrainment 
depends on the near surface stability, which in turn is 
controlled by local atmospheric forcing rf the upper ocean 
layer which nahes a coupled prediction approach necessary. 
A circularity in cause and effect arises because surface 
buoyancy flux induced enrrainment nor only increases the 
depth of the MABL, but it also changes its effect on the 
ocean mixed layer. The role of rhe surface buoyancy flux 
means changes in the oceanic well-mixed layer could be 
either the cause or the result of entrainment. Anorher 
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cause and effect relationship is the forma- ion and dissipa- 
tion of clouds within the ABL. Clouds can be caused by 
changes of the ocean surface temperature. In turn, clouds 
have a profound effect on the short and Long wave radiation 
budget for the OBL. 

These cause and effect relationships show that the 
assessment of relative roles of dynamic processes in borh 
layers requires coupled evaluations using relatively complex 
physical models. A data set collected (Davidson et al 1983) 
over a 48 hour period in May 1978 demonstrates the nature of 
coupled local changes in the atmospheric and oceanic bound- 
ary layers. For a period from 1700 to 0500 local time, 
clouds existed in the HABL and winds increased steadily from 
6 to 9 ms'^ during this period. It is suggested this 
increase was associated with deepening of the MABL (and 

therefore, the entrainment of additional momentum into the 
mixed layer) which is, in turn, associated with the surface 
temperature. The oceanic boundary layer deepened as the 
wind reached a maximum, resulting in the surface temperature 
dropping from 14.5° to 13.0°C. 

There were other atmospheric response and forcing fea- 
tures observed during the entire 48 hour period from 
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5/19/1200 to 5/21/1200 PST which are believed to have been 
associated with coupling between the two layers. Atmos- 
pheric and oceanic mixed layer observations appear in Figure 

4 . 

These additional features are; 

a) The ABL depth increased from 250 to 750 meters over 
the 48 hour period. The changes occured in relatively 
short intervals, from 20/0000 to 20/0400, and from 
21/0000 to 21/0500. The ABL depth remained at near 
constant depths during the intervening 12 to 18 hour 
periods. 

b) The surface layer temperature increase from 20/0000 
to 20/1500 is indicative of entrainment of overlying 
warm air. This entrainmen-^ presumably was a factor 
in the wind speed increase over the same period. The 
speed increase, in turn, presumably caused the warm 
shallow OBL to be destroyed which dropped the surface 
temperature after 20/1800. 

c) The most dramatic change in the MABL occured during 
the period from 21/0000 to 21/0500 which was immedi- 
ately after the surface temperature had dropped due 
to wind mixing. 
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Figure 4. Atmosoher ic and Oceanic Mixed Layer Observations 
durind CEWCOM-78. (a) porential temDerature 
profiles ( t) accusric sounder record' (c) 10 merer 

wind speed (d) 10 meter and surface temoerature 

and (e) XBT rraces. 
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Questions which could be posed on the basis of the above 



data are: 

a) Did the dramatic change in rhe ABL depth at 21/0000 
occur as a result of the ocean forcing the atmosphere? 

b) Does the fact that the mixed layer evolved with rapid 
changes followed by equilibrium periods indicate feed- 
back between the boundaries? 

c) What wculd have occurred if the wind had reached its 
maximum value at a time other than after sunset? 

d) Was the time of occurrence of the wind maximum and, 
hence, the ocean mixing conrrolled by the coupled 
mixed layer system? 

This thesis will attempt to answer some of these ques- 
tions on the basis of the formulated coupled model. As 
seated earlier, all answers can only be obtained through- 
considerable interpretive efforts. 

The scope of the overall scientific effort is illus- 
trated in Figure 5. The local oceanic descriptions will 
made by a bulk oceanic mixed layer model (the Garwood 
Model). The atmospheric descriptions will be mads by the 
Naval Postgraduate School (NPS) Marine Atmospheric Boundary 
Layer model (Davidson et al 1983). Tactically significant 
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operational descriptions (forecasts) to which the improved 
understandings would apply, appear on the extreme right hand 
side of Figure 5. 

B. MARINE ATMOSPHERIC BOUNDARY LAYER (MABL) MODEL 

The NPS boundary layer modal is a zero-order, two layer, 
integrated mixed layer model. The two layers consist of the 
well-mixed, turbulent boundary layer, and the relatively 
non-turb ulenr free atmosphara. The model is based on 
entrainment energetics formulated by Stage and Businger 
(1981) and on radiative transfers described by Davidson e- 
al (1983). 

The inversion (or transition zone) separates one layer 
from the other. In a zero-order model, the inversion is 
assumed to have zero thickness, hence, a discontinuity, or 
jump occurs in the profiles of the ccnsarvative parameters 
at the inversion. 

The present model requires as input, an initial atmos- 
pheric sounding, the mean winds at a level within the sur- 
face layer (10-30 meters) and the surface temperature. 
Boundary layer wind values for the duration of the forecast 
period can be obtained from larger scale models. In this 
study, actual hourly wind speed measurements were input as 
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part of the model initialization. Sea-surface temperature 
remains unchanged from the initial input value. However, in 
the coupled model, the sea-surface temperature will change 
during the forecast period and will be that predicted by the 
OBL model. Well-mixed temperature and humidity are predic- 
ted sc the surface wird and wind shear at the inversion are 
the only atmospheric variables which have to be prescribed. 
The larger scale subsidence, obtained from the large scale 
model must also be prescribed for the forecast period. 

Three methods can be used to compute the large scale 
vertical velocity (subsidence) from single station observa- 
tions. These three methods are well described by Gleason 
(1982) and include the kinematic and adiabatic raerhods and 
-he integration of the moisture budget equation (Q-methcd) . 
Gleason's study shewed the Q-method displayed the most merir 
as a single-station assessment of subsidence. This method 
was used to compute subsidence for this study. The solar 
zenith angle used in the computation of the short wave radi- 
ative flux is determined from the latitude, julian day and 
time of day of the initial inpur. 

The model predicts the time evolurion of the inversion 
height, and the mixed layer values of 9 (temperature) and g 
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(humidity), as well as their jumps, at 33-minute inrervals. 
The prediction model is run out to 24 hours. The steps in 

the prediction computation are shown in Figure 6, where it 
is noted tha-^ procedures are the same for clear and cloudy 
cases except for entrainment computation and estimation of 
cloud top cooling. 

The model uses the standard integrated rate equations 
(Tennekes and Dreidonks, 1981) to predict the time evolution 
of the conservative variables and their jumps at the inver- 
sion. These equations are: 



h(Dx/Dt)= (w'x')o “(w’x*),^'*' source (1) 

h(D^x/Dt)= hrx(ah/dt) - (w'x')q+ (w'xMh* source (2) 



wher 
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Stage and Businger's (1981) entrainment velocity parame- 
rerizacion is used to close the system of equations and 
determine the time evolution of the inversion height. The 
closure assumption states that the dissipation rats of the 
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Figure 6. Schematic cf Input, Prescription and Computing 
Steps in «ABL Prediction, 
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turbulent kinetic energy (TKE) is a frar’rion (1-A) of the 
preduction rate. ' A • is the entrainment coefficient, taken 
as 0.2. 

Mathematical modeling of radiative flux transfer is an 
extremely complicated but still developed subject. Unfortu- 
nately, even the simple radiative transfer models are srill 
extremely complex in comparison to most other physical 
parameterizations used in the mixed layer model. Uncer- 
tainty in background aerosols, atmospheric absorbing gases 
(water vapor, carbon dioxide, nzone) and cloud droplet size 
spectra are sources of error in radiation calculations. 
Long and short wave radiation fluxes are computed sepa- 
rately. 

The long wave radiation flux calculari on was modified to 
permit non-black stratus clouds by introducing the cloud 
emissivioy, £ , which is a function of the total cloud liq- 
uid content, W. Since the cLoud liquid waoer profiles are 
approximately linear with height (Davidson et al, 1982), 
cloud liquid water contenr and cloud emissivioy are given 
by : 



W= 0.5j3a(h - Z^) q,^ 


(3) 


c = 1 - exp (-aW) 


(4) 
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where J?a the density of air ( 1.25X 10’^ gm cm‘^ ) , h is the 
height of the mixed layer (cloud top). Sc is the lifting 
condensation level (cloud bottom), a=0.153 (Slingo, 1982), 
and is the liquid water content at the cloud top. 

The long wave cloud top net radiation flux, can be 

calculated from rhe cloud top temperature, , and the 

effective radiative sky temperature, Tg^y , using the Stefan- 
Boltzmann Law. At the cloud bottom, this flux, Lp^, is cal- 
culated in the same manner using the sea surface 

temperature, Tg , and cloud bottom temperature, Tc . These 
fluxes are given by: 



inh = ^:cCr(V- Tgk/) (5) 

Lnc = £c<r(Tg'' - ) (6) 

where CT is Stefan's constant (h.61 X ) and £1^, , cloud 

emissivity, is obtained from equation 4. The net long wave 
flux at the surface, ?|ong » becomes: 



J’long = ^c) ^sky'l (7) 

where T is the average temperature of the cloud. For the 
cloud free case, the net fluxes are calculated at Z=h and 
Z=0 using the standard method (Fleagle and Businger, 1978) 
of integrating the flux emissivity profile. The net long 
wave flux at the surface for the clear sky case is: 
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P = F 

^long 



( 8 ) 



- F. 



where ?,j is the upward radiative flux and F^j is the downward 
radiative flux. 

Surface flux para me terizario n (bulk aerodynamical) for- 
mulae 



U*= Cj U^o 

T*= cj^^(9o- 9) 
1/2 

q*= (go" 



(9a) 

(9b) 

(9c) 



are used to determine the surface fluxes of momentum, sensi- 
ble heat and latent heat. Th=se fluxes are given by: 



u • w' = u ♦ 

T * w * = u ^ ^ 
q I w' = u *q* 



(m omen turn) 
(sensible heat) 
( 1 at an t heat) 



(10a) 

(10b) 

(10c) 



where Cq and are rhe ten meter stability dependent drag 

coefficient, 9 is zhe potanzial zemperature, and g is zhe 
specific humidity. The subscript zero denotes the surface 
va 1 ue . 

The shorz wave radiative flux is caLcula-^ed using the 
delta-Edingtcn method (Joseph, 1976) . This was added to the 
model zo acccunz for hearing of the mixed layer by solar 
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radiation. Short wave extinction is dominated by scattering 
(as opposed tc absorption). This scattered radiation forms 
a second short wave radiative component usually referred to 
as ’diffuse’ solar radiation. Scattering in the mixed layer 
is due to atmospheric particles, cloud droplets in the 
cloudy case and sea-salt aerosols in the cloud-free case. 
This MABL model computes both the direct and diffuse radia- 
tive components to determine rbe total short wave radiation 
flux at the surface. An explanation of the delta-Edington 
method and all of the parameters, atmospheric factors, and 
equations involved is a very complex subject and is beyond 
the scope of this paper. An excellent: review of this sub- 
ject has been published by Fairall (1981). An important 
short wave radiative parameter that, is prescribed in the A3L 
model is the fraction of reflected short wave radiation, Ag, 
from the sea surface. It is taken as 0.1. 

C. OCEANIC BCONDARY LAYER (OBL) MODEL 

Garwood (1977) developed an ocean mixed layer model 
using the Na vier-Stokes equation of motion with the geostro- 
phic component eliminated, the continuity equation in incom- 
pressible water, the heat equation from the first law of 
thermodynamics, the conservation of salt equation, and a 
linearized equation of state. 
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Prediction of the rate of deepening (ar retreat) of the 
mixed layer is dependent upon an understanding of the dynam- 
ics of the entrainment process. The turbulence of the over- 
lying mixed layer provides the energy needed to destabilize 
and erode the underlying stable water mass (Garwood 1977). 
Therefore, the turbulent kinetic energy budget is the basis 
for the entrainment. This system is closed using a mean- 
turbulent-field modeling of the verrically integrated equa- 
tions for the individual TKE components, plus the inclusion 
of the bulk buoyancy and momentum equations. 

Separate vertical and horizontal equations for TKE are 
used to better handle the mixing process. Buoyancy produced 
energy is somewhat more efficient than shear production as a 
source of energy for vertical mixing because of its unique 
effect on the vertical component of che turbulent velocity. 
The buoyancy equation is generated from the heat and salt 
equations together with an equation of state: 




( 11 ) 



and buoyancy is given by: 



b= q{ fo - f) /fo 



( 12 ) 



32 



where 0 is temperature, S is salinity, j? lensity, g gravity, 
ar.d the constants oc and ^ are rhe expansion coefficients for 
heat and salr, respectively. The -^ilde represents the total 
instantaneous value and the subscript zero denotes an arbi- 
trary, but representative, constant value. Using b instead 
of 9 (temperature) for the definition of buoyancy allows 
this model to be applicable in situations where evaporation 
and precipitation contribute significantly to the surface 
buoyancy flux. Temperature is in most cases considered to 
be nhe dominating factor in the density profile. The rela- 
-ive importance of the salinity effect on the short term 
evolution of the density profile is generally not signifi- 
cant (Miller, 1976). 

Temperature and salinity profiles, and fcr extended 
forecasts, the wind driven horizontal current profiles, are 
required as input for model initialization. These are used 
to compute the mixed layer depth. The Sarwocd OBL model 
defines the mixed layer depth, h , as the shallowest depth 
at which the observed density value, (T\ , is 0.02(T[ units 
greater than the observed surface density value. 

Ocean environmental parameters that have to be pre- 
scribed include; the fraction of short wave radiation 
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absorbed in the top one meter of the ocean. 



the radiation 



extinction coefficient for determination of the downward 
turbulent heat flux, and the critical Richardson number for 
a stability adjustment at the bottom of the mixed layer. 

The parameters required for surface boundary condition 
computations include wind speed and direction, cloud cover, 
sea surface temperature, air temperature (dry bulb), dew- 
point temperature, incident solar radiation, and the rate of 
evaporation (E) and precipitation (P) . 

The turbulent fluxes of latent heat, Qg, and sensible 
heat, , are estimated using the bulk aerodynamic formu- 
lae : 

Qe= Cd (.98 Es - '^io 

Qh= Cd (Ts - T3 )u^o (13b) 

and the net back radiation is esnimated from the empirical 
equation (Husby and Seckel, 1975): 

Q^= 1.14X 10'^(273. 16+Tg)'* (.39 - .5E3^^)(1 - .6C^) (13c) 

where Eg is the sarurated vapor pressure of the marine air 
(0.98 corrects for salt defects) , is the vapor pressure 

of air based on the dew-point temperature, air tempera- 

ture, Tg sea surface temperature, and C is the fracnicnal 
cloud cover. The upward heat flux, Q^, is then given as: 
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Q = Q + 2^. 

*u '^b 



( 14 ) 



The solar radiation, Qg, is givsn by: 

Qg= (1 - aoc*") (1 - .66C^) Qq (15) 

where Qq is the clear sky radiation given by (Seckel and 
Beauday, 1 973) : 

Qq= Aq+ CCS sin ^ + Agcos 2^ + BgSin 2^ (16) 

The constants a and b are adopted from Tabata (1964) and the 
cubic cloud cover correction from Laevas.u (1960). The 

variable is the mid-day elevation angle of the sun. The 

coefficients (Aq,A^, ect.) were calculated by harmonic rep- 
resentation of the values predicted in the Smithscnian 
Ket eorological Tables (list 1953) with 

(2TT/365) (t- 21) (17) 

where t is the julian day of -he year. 

Not all of the incoming short wave radiation penetrates 
rhe ocean mixed layer. Approximately half (for open ocean) 
is absorbed within the first meter. The amount absorbed 
varies from region to region depending on the amount of 
absorbing particles such as phytoplankton and yellow sub- 
stance. More radiation will be absorbed in coastal regions 
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than in the open ocean regions. This portion of absorbed 
short wave radiation is therefore considered to be part of 
the upward heat flux. The ramaindec of the short wave radi- 
ation does penetrate the mixed layer and is attenuated in an 
exponential fashion depending on water turbidi-^y. If the 
fraction of solar radiation absorbed in the first meter, Rf, 
is given, then the net heat at the surface is given by: 



Qnet = Qu ^ (Rf)Qs - Qs (13) 

From these, the surface fluxes of buoyancy (heat and 
salt) and momentum can then be computed. The turbulent tem- 
perature, salinity, velocity and buoyanoy fluxes are given 
by : 



respecti vely , 
lated by; 



(T'w*) = 


"Qnet 


(19a) 


(S*w») = 


(P-2) So 


(19b) 


(u' w* ) = 


u*^ 


(19c) 


(b*w*) = 


q[o( (T ' w' ) - (^(S* w' )] 


(19d) 


where u*. 


friction velocity in air. 


is calcu 


"^s “J^a 


2 

^10 


(20) 



u*= (21) 
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The subscripr zero refers to the surface value and ZTs is the 
surface stress (dynes cm ). In this nodel, a positive 
buoyancy results when Qnet<0 and E>P. Daring daytime peri- 
ods, negative buoyancy flux results from the domination of 
the solar heating at the surface, while during the night 
periods, a positive buoyancy flux is produced by the com- 
bined effect of net long wavs radiation and the upward tur- 
bulent fluxes of heat and moisture. 

New ocean temperature, salinity and wind driven current 
profiles are predicted at one hour intervals. These are 
then used to predict the new mixed layer depth. The steps 
in the prediction computation are shown in Figure 7. 

D. FORMULATED COUPLED 30UNDARI LAYER MDDEL 

The physical processes described in the preceding two 
sections were not alt'^red during the coupling steps. Some 
modifications were involved in Transferring parameters 
between rhe models. Also necessary was the uniformiry of 
units. The only parameter supplied from the ocean model for 
use in atmospheric computations is the predicted sea surface 
Temperature. However, this predicTion is influenced by many 
armospheric parameters. 
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Figure 7 . 



Schematic of Input, Prescription and Computing 
Steps in OBL Prediction. t- y 
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Parameters in the atmosphere initialization include; 
wind speed and direction, dew-point temperature, latitude, 
julian day, and time of day. Short wave radiative flux 
(direct and diffuse) , long wave radiative flux, sensible and 
latent heat fluxes, and the surface momentum flux were com- 
puted by the atmospheric section of the coupled model and 
input to the ocean model. 

The first problem associated with the coupling involved 
^he different prediction time steps of the seperate models. 
The atmospheric model uses a 30-minute time interval, 
whereas, the ocean model uses a one-hour time step interval. 
Assuming no significant sea-surface temperature change would 
take place over the short one-hour interval, the ocean mod- 
el’s prediction scheme was inserted into the atmospher® 
model as a subroutine. This prediction scheme is only 
called every other 30-minute time step. The ocean model's 
initialization routines are executed immediately upon the 
completion of the atmosphere model's initialization. 

The second problem was with wind input. The ABL model 
does not require wind direction as an input parameter, only 
~he magnitude of the wind speed. However, the ocean model 
requires wind direction to compute the horizontal components 
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cf the winds. These components are used to determine the 

2 

ocean turbulent velocity flux, ' components used in the 

momentum budget compu taticns. k new subroutine was added in 
the initialization scheme to determine the horizontal wind 
components from wind speed and direction input. These com- 
ponents are determined by: 

’^lOx = ■ *^10 (22a) 

'^I 0 y = -(cos 9) Uio (22b) 

where u^q specified wind speed at 10 meters and the 

east-west and north-south components are ’^lox '^loy ' 

respectively. 0 is the direction from which the winds are 
blowing relative to true north. Use of equations 20 and 21 
and the following relationship: 

■Cs= (23) 

2 

allows the ocean turbulent velocity flux ) cc be dener- 

mined (assuming rhe density of seawater, , to be approxi- 
mated by 1) by first computing the wind stress from 

atmospheric parameters and then equating ‘rhe ocean turbulent 

velocity flux to the computed wind stress value. The hori- 

, 2 
zontai components of the turbulent velocity fluxes u*^^ 

2 

u*^y , are then easily computed by: 
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( 24 ) 



U*^=Cfl)U,o^= C^(U,ox^ + U,oy^) 

u=«Wx^= C(,Uio/ (25a) 

^(J^lOy (25b) 

The atmospheric derived values are computed using the HKS 
unit system and had to be converted to the CGS system prior 
zo any oceanic computations. This conversion is done at the 
start of each ocean prediction run. 

As discussed in section B, the short wave radiative flux 
is calculated using the delta-Sd ingt on method. The long wave 
radiative, s<=‘nsible heat and latent heat fluxes are computed 
using equations 10a, 10b and 10c, respectively. A unit con- 

version from Watts m to cal cm s is required on these 
values prior to use by ocean model. This conversion is 
again done a~ the start of the ocean prediction run. A flow 
diagram of the steps in the coupled model’s prediction com- 
putation is shown in Figure 8. 
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Figure 8 



Scnematic of Atmospheric and Oceanic 
Hoasl. The. ocean model is called at 
oO-mnnure time step to forecast prescribed 
surface -erperarure for the MA3L orediction. 



Coupled 
eyftY othe 



42 




III. DATA AND MODEL RESULTS 



The data sets used in this study are from the Coopera- 
tive Experiment on West Coast Dceanography and Meteorology 
from 21 September to 12 October 1976 (CEWCOM-76) and from 5 
to 23 May 1978 (CEWCOM-78) . Both were conducted in the Los 
Angeles-San Diego basin. Oceanic and atmospheric data were 
collected simultaneously throughout these time periods. 
XBT ' £ were deployed from the R/V ACANIA on a regular basis 
(2-3 per day) in conjunction with atmosphere radiosondes. 
Sea surface temperature and atmospheric data such as winds, 
temperature, and relative humidity were also obtained con- 
tinuously through these periods on board the R/V ACANIA. 

Three 2U-hour periods were selected for model initiali- 
zation and verification. The first two cases cover a 
48-hour period, 3/19/0500 to 5/21/0500 PST, during 
CEWCOM-78. This period coincides with the observarions dis- 
cussed in chapter 2, in which coupled air-sea interactions 
were cited as possible explanations for changes in both the 
atmospheric and oceanic mixed layer. Satellite imagery 
showed uniformly increasing srratus coverage (thin zo heavy) 
during rhis 48-hour period. A clear sky situation was 
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chosen fox the thir5 case. This 24-hcur period, 10/04/1900 
to 10/05/1900 PST, was during CEWCOM-76. 

The general location of the R/V ACANIA during these 
periods is shown in Figures 9, 10 and 11. During case 1 the 
ship was at anchor for the majority of the time period. The 
anchorage location is shown in Figure 9. For cases 2 and 3 
the R/V ACANIA was cruising slowly (2-3 knots) re-^urning to 
an initial point approximately every 12 hours. Winds were 
primarily from the west during the case 1 period. The 
winds shifted slightly to a west-northwest direction midway 
through the case 2 period. Winds for case 3 were primarily 
from the northwest and were much slower than in either case 
1 or 2. 

Two predictions were analyzed for each case, one predic- 
tion was based on the coupled air-sea boandary layer model, 
and the second prediction was with the atmospheric model 
with a fixed surface temperature. A separate oceanic pre- 
diction using fixed atmospheric parameters (persistence) was 
not made. This was the result of an ini'ial modification 
that removed the wind stress and heat flux computation pack- 
age in the ocean model in preparation for coupling with the 
atmosphere model. 
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Figure 9. Anchorage Position of F/V ACANIA in the Vicinity 
of San Sicolas Island froai 5/19/78 to 5/20/78. 

Presentanion of results will be in two parts. The first 
to be presented is the coupled model’s performance with 
respect to observed oceanic and atmospheric variations. The 
second will deal wi~h a comparison of predicted values 
between the coupled model and the atmospheric model. All 
rabies and figures are grouped together by case number at 
the end of the chapter. 
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Figure 10, Positions of R/V acANIA from 5/18/78 to 5/25/78. 
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Positions of R/V ACANIA from 10/4/76 to 10/6/7 
Horizontal tamperatare structure is shown by t 
dashed lines. 



A. RESULTS FROM COUPIED MODEL 

1. CASE 1 (0 500 Max 0 500 20 Ha^) 

A listing of predicted atmospheric and oceanic mixed 
layer quantities for this case appear in Table I. Values 
are for each 30-minute time step and since oceanic values 
are only computed at one-hour intervals, the same value 
appears at two time steps. The parametsrs listed include: 

height of inversion (Zi) , height cf lifting condensation 

level (Zlcl) , wind speed (Wind), atmospheric mixed layer 
potential temperature (Th) , evaporation duct h^^ight (Ze) , 
sea surface temperature (SST) , ocean mixed layer depth 
(MLD) , and the net ocean surface heat flux (TW) . 

Predicted values are also shoi^n in Figures 12 
through 17. Figure 12 shows the time evolution cf predic- 
"■ed sea- surface temperatures (solid line) and the values of 
the measured surface temperatures (dashed line). Figure 13 
shows the time evolution of the ocean mixed layer depth. 
Figures 14-16 show ocean vertical profiles for temperature, 
u and V wind driven current velocities, and the salinity 
change from an ini'^ial value of 33.5 ppt, respectively. 
Figure 17 shows predictions of atmospheric parameters. 
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High solar insulation luring tha daytime rasulrs in 
a net downward heat flux (ocaan warming) . This can be read- 
ily observed in the net surface heat flux values (TW) , in 
Table I, with the negative sign appearing during the daytime 
hours. Maximum heating of the ocean occurred between 1300 
and 1400 local time. This same pattern, warming in the day- 
time hours, is also observed in the sea surface temperature 
values. Table I and Figure 12. 

The winds appeared to be fairly constant, 8ms*’, 
between 0500 and 1100 on the nineteenth. Sind mixing coupled 
with convective cooling at the surface caused the mixed 
layer depth to deepen slightly from 19.5 to 20.9 m (Figure 
13) . This deepening and net heat loss at tha surface caused 
the surface temperature to drop. As heating of the sea sur- 
face began to increase due to the increase of solar radia- 
tion, a shallow mixed layer of 1 4 meters formed. 

The winds rapidly decreased from 3 to 5 ms*’ at 1200. 
This coincided with a near maximum value of radiative hear- 
ing causing a shallow 7 merer mixed layer ~o form in less 
■chan 1 hour. As daytime heating continued, a shallower 
mixed layer of 3.5 meters formed. This process is referred 
to as the ’afternoon effect’ in acoustic Literature and is a 



49 



quite common phenomenon in many areas of the ocean. In Fig- 
ure 14, the 7 meter mixed layer is barely discernible in the 
1300 vertical temperature profile. The weak gradient at the 
base of the mixed layer only required weak mixed layer ero- 
sion (increase in winds or convective cooling) to break down 
the shallow mixed layer. During the late afternoon, the 
winds began to increase, and with the decreasing of solar 
radiation, the shallow thermocline was eroded away. The MLD 
continued deepening for the remainder of the time period due 
to wind- driven and convective mixing. 

Sines sea-surface temperature is the sole input from 
the ocean model to the atmosphere model, a comparison is 
made of predicted and actual sea-surface temperature values 
(Figure 12). Since the R/7 A3ANTA was at anchor during the 
majority of this time period, this is considered the best 
case to determine how well sea surface temperature changes 
due to atmospheric forcing are being predicted. The coupled 
model predic"*ion yielded a 0.25°C temperature increase dur- 
ing the daytime heating period with a maximum temperature of 
15.4^0 at 1600 local time. A maximum surface temperature of 
15.6°C was observed at 1600 givina a 0.2*^C error. 
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After the daytime heating period, the surface tem- 
peratures began to drop. The minimum predicted surface tem- 
perature at the end cf the 24-hour period was 15.2°C which 
yields a 0. 2°C temperature drop. The 0.06°C temperature 
drop prior to ocean warming makes the magnitude of the pre- 
dicted diurnal warming and cooling to be approximately 
equal. The actual minimum temperature at the end of the 
period was 14.96°C. Although there are small errors in the 
magnitude of the diurnal heating and cooling affects, the 
time evolution of this phenomenon is predicted quite well. 

Precipitation was not observed during this period, 
therefore the surface moisture flux and surface saliniry 
flux were due entirely to evaporation. The slight predicted 
increase in salinity from 33.5 to 33.5038 ppt, shown in Fig- 
ure 16, is due to evaporation. 

An unstable condition existed in the atmosphere dur- 
ing this prediction period. A strong bnoyancy flux at the 
surface increases entrainment and should cause the inversion 
height to increase unless subsidence (associated with shal- 
lowing of the inversion height) is strong enough to prevent 
this increase. In this case, subsidence was relatively weak 
(-. 0028 ms"^) and entrainment was able to increase the inver- 
sion height. 
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Cloud occuran ce/dissipat ion was of interest because 
of its effect on radiation at ttie surface. With the lifting 
condensation level, LCL, lower than the inversion throughout 
this time period, stratus clouds should be present during 
rhe entire period. This feature was observed. During this 
period the LCL decreased from 243 to 125 meters while the 
inversion height increased from 256 to 348 meters. This 
would indicate a thickening of the clouds over this 24-hour 
period which was also observed to be true. The air-sea 
coupling effect on predicting these cloud features cannot 
be identified until the results with a constan"^ sea surface 
-emperature forcing are examined. However, on the basis of 
the predicted temperature change, a significant influence by 
the coupling is not expected. 

2. CASE 2 (0 500 20 May to 0 500 21 May) 

A listing of predicted atmospheric and oceanic mixed 
layer quantities for this period are given in Table III and 
shown graphically in Figures 19 through 24. As in case 1, 
daytime heating at the surface is apparent by the negative 
signs in front of the net surface heat flux values. The 
maximum downward heat flux value was not as high as in case 
1 and was probably due to thicker cloud conditions for this 
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case. However, increasing surface temperatures were 
observed during this heating period. 

The winds increased steadily from 6 to 11 ms’’ 
throughout the first half of this period and slowly 
decreased to about 9.4 ms*’ at the end of the period. The 
higher winds should be associated with increased wind mixing 
of the ocean mixed layer. The increased winds caused the 
mixed layer depth to deepen at the end of this peiod. This 
also reduced the shallowing of the mixed layer during the 
afternoon radiative heating period. The shallowest mixed 
layer formed at 16 meters versus 3, 5 meters in the previous 
case. In Figure 21, the 20 meter thermocline is seen to be 
very weak, and with the higher winds, the shallower after- 
noon mixed layer lasted only a couple of hours. 

The effects of the diurnal hearing and cooling peri- 
ods caused a 0.1°C temperature increase and a 0.15°C temper- 
ature decrease, respectively. Borh values are also smaller 
than those predicted in case 1 . Beasons for the smaller 
predicted increase are; 1) thicker clouds reduced the amount 
of solar hea-^ing during the day, therby reducing the diurnal 
xemperature increase; 2) the higher winds increased the ver- 
tical mixing, causing cooler water below to be mixed into 
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the upper layer and also increased the latent heat loss due 
to evaporation, thus reducing the total net heat flux at the 
surface; and, 3) the larger air-sea temperature difference 
increased the sensible heat loss, which also reduced the net 
heat flux. The decrease in the downward solar radiative 
heating and increase in sensible and latent heat losses were 
somewhar balanced by the decrease in the long wave heat loss 
due to the presence of the clouds. This decrease in the 
long wave heat loss was also the primary reason for the 
slightly smaller value in the nighttime cooling period. 

Comparisons cf the predicted surface temperature 
values with the actual values (Figure 19), show a noticeable 
difference. This is probably due to rhe fact than at 0700 
local time on 20 Hay, the R/7 ACANIA began cruising slowly 
to the northwest. The regional horizontal temperature 
sxructure must be taken into account. The southern Califor- 
nia coastal barhymetry and associated currents greatly 
influence the sea surface temperature sxructure. The Cali- 
fornia Current brings in colder water from the north as it 
heads south along the coast. Near Point Conception, the 
current begins to turn to the southwest away from the ccasx. 
Due to coastal bathymetry. Point Conception also provides 
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sheltering for the intrusion of warm water along the coast 
from the south by the Davidson Current. The ‘=ddies and 
meanders of cold and warm warer make the sea surface temper- 
ature structure complex in the CEWCDM experiment areas. The 
surface temperature variation observed by the R/V ACANIA 
appears to be due to oceanographic features. 

The 1®C temperature increase in the first two hours, 
when the ship began moving to the northwest, indicates that 
the ship's "^ransit was through one of the small warm eddy 
features. Two reasons for this are 1) Radiative heating of 
the ocean did not even begin until after the observed sur- 
face temperature had already begun to rapidly decrease, and 
2) the higher winds would tend to mix the surface waters 
with the cooler water from below, thus causing a decrease in 
temperature not an increase. The short time duration of 
this warm feature also suggests a transit through an ano- 
maly. With the ship only cruising at 2-3 knots, the warm 
anomaly still only lasted less than 4 hours. 

The ship turned on a northerly course at 1300 and 
the rapid temperature decrease stopped. It is reasonably 
safe to infer from the sudden change in the surface tempera- 
ture trend (rapid decrease to slight increase) that the ship 
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is initially crossing the temperatare gradient contours and 
then followed a parallel course with rhese contours. This 
would make oceanic conditions more compatible with those in 
case 1. Since the 0.4®C temperature increase during the 
daytime heating period compares favorably with the value 
observed while the ship was at anchor, the increase is 
therefore presumably due to the diurnal heating. 

The ship turned back to a west-northwest course 
shortly after 1800 and the surface temperature immediately 
began to decrease 1.5°C in 4 hours. Two hours after the end 
of this period, 0700, the ship came about and paralleled its 
previous course but in the opposite direction. An examina- 
tion of the observed sea surface temperatures shows a rapid 
increase comparable to the decrease while heading in the 
opposite direction. This further substantiates the observed 
large temperature changes are not due to atmospheric forc- 
ing, but rather to the ship's movement through an oceano- 
graphic eddy or front. 

Unstable atmospheric conditions existed during this 
period as well. The 1.5°C drop in the air temperature and a 
much smaller drop in the ocean temperature (0.2°C) coupled 
with higher wind speeds resulted in a large surface buoyancy 
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flux. This resulted in a much higher predicted inversion 
height change due to entrainment. The slightly smaller sub- 
sidence rate was also a factor in the greater predicted 
increase of the inversion height. 

The lifting condensation level fluctuates up and 
down unlike the steady decrease in height in case 1. How- 
ever, the height of the LCL never increased above the ini- 
tial height and so with the large increase in height of the 
inversion, thickening of the stratus cloud deck is predicted 
to continue during this period as well. 

3. CASE 3 (JliOO 04 OCT to 1900 05 OCT) 

A listing of predicted atmospheric and oceanic mixed 
layer quantities are given in Table V and shown graphically 
in Figures 26 through 31. Although the time period for this 
case is in October, there are still many features occuring 
that are similar to the first two cases in May. The primary 
differences in this case, other than being a cloud-free 
case, are 1) the air and surface temperature are much warmer 
and 2) a much higher value of subsidence was observed. 

Daytime heating is apparent in predictions for the 
clear sky case, as expected. The magnitude of the net ocean 
surface heat flux is, however, smaller than the value for 
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the thin stratus deck condition in case 1. This is presum- 
ably due to the difference in seasonal heating rate and will 
not be considered part of the local coupling effects. 

The winds in case 3 are much weaker than in the 
first two cases. with the light winds and high downward 
surface heat flux during the day, a much shallower mixed 
layer might be expected. The ocean mixed layer depth (MLD) 
values in Table V (also shown graphically in Figure 25) show 
that an ocean mixed layer of 1 meter suddenly forms in less 
•^han an hour after heating began at 0910. This shallow 
thermocline is much stronger than in either case 1 or 2, and 
persists throughout most of the day (Figure 28). As the 
winds slowly increased late in the afternoon, the shallow 
mixed layer deepens. Without winds to mix this shallow 

region, surface heating is much greater in this case. A 
0.6°C surface temperature increase was predicted to occur. 

As in case 2, the ship’s movement played a major 
role in the large variations in observed surface temperature 
values. This prevented the verification of the predicted 
surface temperature increase during the daytime heating 
period. Figure 11 shows the complex horizontal temperature 
structure in the CEWCOH-76 experiment area for the first 
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week of October. From this figure, it is hard to explain 
the large 2.8°C temperature drop during the first half of 
the period as the ship appeared to be heading into warmer 
water. The rapid temperature iacrease, starting at 0900, is 
obviously due to the ship transiting into the warm eddy fea- 
ture north of Santa Catalina Island {Figure 11). 

Unstable atmospheric conditions are also present 
during this time period. The high subsidence caused the 
inversion to be forced downward to within 72 meters of the 
surface. The lifting condensation level decreases even more 
rapidly and finally becomes lower than the inversion at the 
last hour of the forecast period. This implies the fcrma- 
■^ion of a low straxus cloud deck. From Figure 33, fog and 
stratus clouds were observed to form in the CEWCOM experi- 
ment area. 

E. COMPARISON OF COUPLED AND ATMOSPHERIC MODEL 

The second and final phase of this study is ro compare 
the previous coupled air-sea boundary layer model resulxs 
with the atmospheric model results using a consxant sea sur- 
face temperature. This comparison is tn derermine if any 
differences existed between the atmospheric predictions, and 
if so, determine if these differences are due no coupled 
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interactions. A listing of the atmospheric model predic- 
tions of atmosphere mixed layer parameters for cases 1, 2 

and 3 are found in Tables II, 17 and 71, respectively. 
Graphical presentations of these results are shown in Fig- 
ures 18, 25 and 32. All tables and figures are grouped 

together by case number at the and of this chapter. 

The initial comparison reveals no significant difference 
in the two predicted sets of mixed layer parameters, indi- 
cating a negligible effect of coupled interactions on atmos- 
pheric predictions for xhis short time period. The 
predicted inversion heights are essentially identical in 
both sets of predictions for all three cases (maximum dif- 
ference of 2 meters) . This is also foand to be true for 
humidity and even the evaporation duct height predictions. 
This observation is somewhat surprising since there are sig- 
nificant interactions between the surface fluxes (evapora- 
tion) and these parameters. The close proximity of these 
parameters to the ocean surface would also imply a more sig- 
nificant impact even with small changes in sea-surface tem- 
peratures. 

The only parameter that seems to be clearly affected by 
the coupling is the lifting condensation level. However, 
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even for the cloudy conditions of cases 1 and 2, this shows 
only a slight deviation. k maximum difference of only 5 
meters is observed between the predictions for these cases. 
In the clear sky case under warmer and more humid condi- 
tions/ there is a 30 meter difference observed in the two 
LCL predictions. However this occurred as the atmosphere 
approach near neutral conditions. This may indicate the 
model's inability to make reliable forecasts under stable 
atmospheric conditions. 

Since a comparison of predicted ocean parameters between 
a coupled and independent ocean model was not made, no con- 
clusions on this topic can be made at this time. Although 
tentative results indicate no definite advantages using a 
coupled model for the atmospheric predictions, the reverse 
may very well be true for oceanic predictions . 
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IV. CONCLUSION S/R ECOMMENDAnONS 



On the basis of a limited study of several 24-hour pre- 
dictions, surface temperature variatioas are not large 
enough to require ocean mixed layer processes in short term 
atmospheric forecasts. CareCully selected data sets and 
considerable interpretive efforts are essential if general 
conclusions are to be drawn on the relative roles of dynamic 
processes in bo-^h boundary layers. Further studies are 
required to determine the model's ability to better handle 
forecasts at longer time scales. Further studies should 
also include sensitivity analyses to define situations for 
which air-sea interaction effects are greatest. Studies of 
oceanic sensitivity to atmospneric forcing are required to 
determine the impact it will have on acoustical forecasts. 

The CEWCOM-76 and 78 data, although able to supply 
simultaneous and oceanic data, suffered in two major 
respects First, the complexity of the horizontal temperature 
structure with meanders of warn and cold eddy features had 
definite impacts on observed data. Second, data was col- 
lected while the R/V ACANIA was transiring through a complex 
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ocean temperature structure. Changes in observed conditions 
under these circumst ances can and did mask the interaction 
processes of interest. 

Although the primary results of this study indicate 
coupled air-sea models are probably not necessary for short 
term atmospheric forecasts, a coupled model still has advan- 
tages, especially if long term forecasts are to become pos- 
sible. Such a model would not be the single one-dimensional 
one considered here. However, this model offers many poten- 
tial benefits as a multi-purpose model with single station 
assessment and prediction capabilities in the absence of 
shore-based support. A sound veloci-y profile (SVP) predic- 
tion capability could be easily incorporated into this 
model. In conjunction with present prediction capabilities, 
a multi-purpose model would be available for a wide range of 
needs. Ocean mixed layer depth and S7? predictions can be 
used in predicting ASW sensor capabilities. Environmental 
effects on the performance of weapons and sensors that 
depend on electromagnetic and electro-optical wave propaga- 
tion can be evaluated using the same modal. 
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TABLE I 



Coupled Model Predicted Atmospheric and Oceanic Mixed Layer 

Quantities- Case 1. 



HOUR 


Zi 


Zlcl 


WIND 


Th 


SST 


MLD 


TU 


Zs 


0430 


256 


243 


0 . 0 


13.3 


15.2 


0 . 0 


0.0000 


0 . 0 


0500 


252 


216 


8.0 


13.4 


15.2 


19.5 


0 .0000 


7.2 


0530 


250 


174 


8.0 


13.3 


15.2 


20 . 1 


. 0067 


6.9 


0600 


250 


130 


8.0 


13.0 


15.2 


20 . 1 


. 0067 


6.4 


0630 


252 


116 


8.0 


12.9 


15.1 


20 . 4 


. 0 028 


6 . 0 


0700 


256 


121 


8.1 


12.8 


15. 1 


20.4 


.0028 


6 . 0 


0730 


259 


127 


8. 1 


12.9 


15.1 


20 . 7 


.0011 


6 . 1 


0800 


262 


132 


8.1 


12.9 


15.1 


20.7 


.0011 


6 . 1 


0830 


265 


136 


8.1 


12.9 


15.1 


20 . 8 


-.0007 


6.2 


0900 


267 


139 


8 . 1 


12.9 


15 . 1 


20.8 


- .0007 


6.2 


0930 


270 


142 


8.1 


12.9 


15.2 


20 . 9 


-.0028 


6.3 


iOOO 


272 


145 


8.2 


13 . 0 


15.2 


20.9 


- .0028 


6.3 


1030 


274 


148 


8.2 


13.0 


15.2 


16 . 4 


-. 0045 


6.3 


i iOO 


276 


151 


8.2 


13.0 


15.2 


16.4 


- .0045 


6.3 


U30 


278 


153 


6.6 


13.1 


15.2 


14 . 1 


-.0057 


6 . 3 


1200 


279 


153 


5.0 


13.0 


15.2 


14 . 1 


-.0057 


5.7 


1230 


280 


153 


5.1 


13. 0 


15.3 


7.2 


-. 0066 


5.2 


1300 


281 


153 


5.3 


12.9 


15.3 


7.2 


-.0066 


5.2 


1330 


282 


153 


5.4 


12.8 


15.3 


3.6 


-.0066 


5.2 


1400 


283 


153 


5 . 6 


12.8 


15.3 


3.6 


- . 0066 


5.3 


1430 


284 


152 


5.8 


12.8 


15.4 


5. 1 


-.0052 


5.4 


1500 


285 


150 


5.9 


12.7 


15.4 


5 . 1 


-.0052 


5.4 


1530 


286 


143 


6 . i 


12.7 


15.4 


6.5 


-. 0032 


5.5 


1600 


288 


145 


6.3 


12.6 


15.4 


6.5 


-.0032 


5.6 


1630 


289 


142 


7.0 


12.6 


15.4 


8.3 


-.0008 


5 . 6 


170 0 


OO V 

C . / J . 


138 


7.8 


12.5 


15.4 


8.3 


- .0008 


5.9 


1730 


293 


134 


7.7 


12.5 


15.3 


12 . 0 


. 0014 


6 . 1 


1800 


296 


130 


7.7 


12.5 


15.3 


12.0 


.0014 


6 . 1 


1830 


299 


127 


7.6 


12.4 


15.3 


14.8 


. 0 029 


6 . 1 


1900 


301 


125 


7.6 


12.4 


15.3 


14.8 


.0 029 


6 . 1 


1930 


304 


125 


7.5 


12.4 


15.3 


17 . 0 


.0034 


6 . 1 


2000 


307 


125 


7.9 


12 . 4 


15.3 


17 . 0 


. 0034 


6.1 


2030 


310 


127 


8.3 


12.4 


15.3 


18.7 


.0034 


6 . 2 


2100 


313 


128 


8.1 


12.4 


15.3 


18.7 


.0 034 


6.4 


2130 


316 


129 


7.9 


12.3 


15.3 


20 . 6 


. 0 038 


6 . 3 


220 0 


319 


129 


7.7 


12.3 


15. 3 


20.6 


. 0 038 


6.3 


2230 


322 


129 


7.5 


12.3 


15.3 


21 . 1 


. 0 036 


6.2 


230 0 


324 


129 


7.3 


12.3 


15.3 


21 . 1 


.0036 


6.2 


2330 


327 


129 


7.0 


12.3 


15.2 


21 . 3 


. 0 035 


6 . 1 


0000 


329 


129 


6.8 


12.2 


15.2 


21.3 


.0 035 


6 . 0 


0030 


331 


128 


6.5 


12.2 


15.2 


21 . 5 


.0033 


6 . 0 


0100 


333 


127 


6.3 


12. 1 


15.2 


21.5 


.0 033 


5.9 


0130 


335 


127 


6. 1 


12. 1 


15.2 


21 . 6 


. 0 032 


5 . 9 


0200 


337 


126 


6.1 


12.0 


15.2 


21 . 6 


.0032 


5.8 


0230 


339 


125 


6.2 


12. 0 


15.2 


21.7 


.0031 


5.9 


0300 


341 


125 


6.2 


12.0 


15.2 


21.7 


. 0031 


5.9 


0330 


343 


125 


6.2 


11 . 9 


15.2 


21 . 9 


.0 032 


5.9 


0400 


344 


125 


6.2 


11.9 


15.2 


21.9 


.0 032 


6 . 0 


0430 


346 


125 


6.2 


11 . 9 


15.2 


22 . 0 


.0033 


6 . 0 


0 50 0 


348 


125 


6.2 


11.9 


15.2 


22.0 


.0033 


6.0 
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TABLE II 



MABL 


Model 


Predi 


cted Atmospher 
Case 1 


ic Mixed 

• 


Lay er 


Quantit 


ies- 


HOUR 


Zi 


Zlcl 


WIND 


Th 


SST 


MLD 


TU 


Ze 


0430 


256 


243 


0 . 0 


13.3 


15.2 


0 . 0 


0 .0000 


0 .1 


0500 


252 


216 


8.0 


13.4 


15.2 


0 . 0 


0 . 0000 


7.; 


0 530 


250 


174 


8.0 


13.3 


15.2 


0 . 0 


0 .0000 


6 


0600 


250 


130 


8.0 


13. 0 


15.2 


0 . 0 


0 . 0000 


6. 


0630 


252 


117 


8.0 


12.9 


15.2 


0 . 0 


0 .0000 


6 . 1 


07C0 


256 


123 


8.1 


12.9 


15.2 


0 . 0 


0 . 0000 


5 . ' 


0730 


259 


128 


8.1 


12.9 


15.2 


0 . 0 


0.0000 


6 . 1 


0800 


262 


133 


8.1 


12.9 


15.2 


0 . 0 


0.0000 


6. 


0830 


265 


137 


8.1 


12.9 


15.2 


0 .0 


0 .0000 


6.: 


0900 


267 


140 


8.1 


12.9 


15.2 


0 . 0 


0. 0000 


6.; 


0930 


269 


143 


8.1 


13 . 0 


15.2 


0.0 


0.0000 


6 . : 


iOOO 


272 


146 


8.2 


13. 0 


15.2 


0 . 0 


0. 0000 


6 . ; 


1030 


274 


149 


8.2 


13.0 


15.2 


0 . 0 


0.0000 


6 .; 


iiOO 


276 


151 


8.2 


13.0 


15.2 


0 . 0 


0. 0000 


6 . 


1130 


277 


153 


6.6 


13.1 


15.2 


0 . 0 


0.0000 


6 .; 


1200 


279 


153 


5.0 


13.0 


15.2 


0 . 0 


0.0000 


5.; 


1230 


280 


153 


5.2 


13.0 


15.2 


0 . 0 


0.0000 


5.! 


1300 


281 


152 


5.3 


12.9 


15.2 


0 . 0 


0.0000 


5.: 


1330 


282 


152 


5.5 


12.8 


15.2 


0 . 0 


0 .0000 


5.: 


1400 


283 


ISO 


5.6 


12.8 


15.2 


0 . 0 


0.0000 


5. 


1430 


284 


149 


5.8 


12.7 


15.2 


0 . 0 


0 .0000 


5 .' 


1500 


285 


147 


6.0 


12.7 


15.2 


0 . 0 


0.0000 


5 . 


1530 


286 


145 


6 . 1 


12.6 


15.2 


0 . 0 


0 . 0000 


5 .. 


1600 


288 


141 


6.3 


12.6 


15.2 


0 . 0 


0.0000 


5 . ‘ 


1630 


289 


137 


7 . 0 


12.5 


15.2 


0 . 0 


0 .0000 


5 .! 


1700 


291 


133 


7.8 


12.5 


15.2 


0 . 0 


0.0000 


6 . 


1730 


294 


130 


7.7 


12.5 


15.2 


0 . 0 


0.0000 


6.: 


1800 


296 


126 


7.7 


12.4 


15.2 


0. 0 


0.0000 


6. 


1830 


299 


123 


7.6 


12.4 


15.2 


0 . 0 


0.0000 


6 .: 


1900 


302 


122 


7.6 


12.4 


15.2 


0 . 0 


0.0000 


6 . ' 


1930 


305 


122 


7.5 


12.3 


15.2 


0 . 0 


0.0000 


6 . 


2000 


308 


123 


7.9 


12.3 


15.2 


0 . 0 


0.0000 


6 . 


2030 


311 


124 


8.3 


12.3 


15.2 


0 . 0 


0.0000 


6 .: 


2100 


314 


126 


8 . 1 


12.3 


15.2 


0. 0 


0.0000 


6 . 


2130 


317 


126 


7.9 


12.3 


15.2 


0 . 0 


0.0000 


s . ■ 


2200 


320 


127 


7.7 


12.3 


15.2 


0 . 0 


0.0000 


6 . 


2230 


1 o -j 


127 


7.5 


12.3 


15.2 


0 . 0 


0 .0000 


6.: 


2300 


TOC 


128 


7.3 


12.3 


15.2 


0 . 0 


0.0000 


6 . 


2330 


328 


128 


7 . 0 


12.2 


15.2 


0 . 0 


0.0000 


6 . : 


0000 


330 


127 


6 . 8 


12.2 


15.2 


0 . 0 


0.0000 


6 . 


0030 


332 


127 


6 . 6 


12.2 


15.2 


0 . 0 


0.0000 


6 . 


0100 


335 


127 


6 . 3 


12.1 


15.2 


0 . 0 


0.0000 


5 . 


0130 


337 


126 


6 . 1 


12. 1 


15.2 


0 . 0 


0.0000 


5 .’ 


0200 


338 


125 


6 . 1 


12. 0 


15.2 


0. 0 


0.0000 


5. 


0230 


340 


125 


6.2 


12. 0 


15.2 


0 . 0 


0.0000 


5 . 


0300 


342 


125 


6.2 


11.9 


•< c. o 


0 . 0 


0.0000 


5 . 


0330 


344 


125 


6.2 


11 . 9 


15.2 


0 . 0 


0.0000 


6 . 


0400 


346 


125 


6.2 


11.9 


15.2 


0 . 0 


0 . 0000 


6 . 


0430 


347 


125 


6.2 


11 . 9 


15.2 


0 . 0 


0.0000 


6 . 


0500 


349 


125 


6.2 


11.8 


15.2 


0. 0 


0.0000 


6. 
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Figure 12. Time Evolution of the Sea-Surface Temperature- 

Case 1. Predicted values are shown by the solid 
line. Observed values are shown by -he dashed 
line . 




Figure 13. Time Evolution of Ocean Mixed Layer Depth (MLD) - 
Case 1 . 
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Figure 14. Selected. Ocean Vertical Temperature Profiles- 
Case 1 . 



0 



-10r 



X 

■*» 

a 

5 -28 



-30 



Proftlo times 
1380 2030 

Tin] rr IT I rTTj]*rrrr|TnT]YrTyrn'i]TnT 






iiuluu 



LLuiLl 



U 



f 

t! 



400 

mTjTrrjpTTrpTTT 

I 
1 



V. 



-20 -10 0 10 -20 -10 0 10 -20 -10 0 

Current (m/s) 



10 20 



Figure 15. Selected Ocean Vertical Profiles of the 

Horizontal Wind Driven Currents- Case 1. Fast- 
west component u (solid line) North-south 
compcneni: v (dashed) . 
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Figure 16. 



Selected Charge in Dcaan Vertical Salinity 
Profiles- Case 1. 
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Figure 17. Graphical .Disolay of ihs Coupled Mod®l 
Axmospheric Mrxsd Layer Quantities- Cas 
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Figure 18. Graphical .Display of the MA3L Model P 
Armospher zc Mixed Layer Quantities- C 
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TABLE I II 



Coupled Model Predicted Atmospheric and Oceanic Mixed Laye: 

Quantities- Case 2. 



HOUR 


Zi 


Zlcl 


WIND 


Th 


SST 


MLD 


TU 


Ze 


0430 


351 


316 


0 . 0 


12 . 0 


15.0 


0 . 0 


0 . 0000 


0 . 0 


0500 


350 


284 


6.4 


12. 0 


15.0 


22 . 5 


0.0000 


7.6 


0530 


351 


249 


6.4 


11.8 


14.9 


22.7 


.0 057 


7.3 


0600 


353 


222 


6.4 


11.7 


14.9 


22.7 


. 0 057 


7 . 0 


0630 


357 


208 


6.4 


11.7 


14.9 


22.8 


.0032 


6.8 


0700 


362 


205 


6.5 


11.7 


14.9 


22.8 


. 0032 


6.8 


0730 


368 


209 


6.7 


11.7 


14.9 


22 8 


.0019 


6 . 8 


0300 


373 


215 


6.9 


11.8 


14.9 


22.8 


.0019 


6.9 


0830 


378 


222 


7.2 


11.9 


14.9 


22 9 


.0004 


7 . 0 


0900 


383 


229 


7.5 


12.0 


14.9 


22 .’9 


.0004 


7.1 


0930 


388 


235 


7.8 


12.1 


14.9 


22.9 


- .0014 


7.2 


1000 


393 


242 


8.7 


12.2 


14.9 


22.9 


-.0014 


7.4 


1030 


399 


243 


8.5 


12.3 


15.0 


20 . 1 


-.0028 


7.7 


1100 


403 


254 


8.3 


12 . 4 


15.0 


20 . 1 


-. 0028 


7.7 


1130 


408 


260 


8.2 


12.5 


15. 0 


20 . 0 


- .0038 


7.6 


1200 


413 


264 


8.1 


12.6 


15.0 


20. 0 


-. 0033 


7.6 


1230 


417 


268 


8.1 


12.7 


15. 0 


16.8 


- .0044 


7.6 


1300 


421 


272 


8.2 


12.7 


15.0 


16.8 


- . 0044 


7.5 


1330 


425 


275 


8.3 


12.8 


15.0 


17.3 


- .0042 


7.6 


1400 


430 


277 


8.5 


12.3 


15.0 


17.3 


- .0042 


7.6 



1430 

1500 

1530 

1600 

1630 

1700 

1730 

1800 

1330 

1900 

1930 

2000 

2030 

2100 

2130 

2200 

2230 

2300 

2330 

0000 

0030 

0100 

0130 

0200 

0230 

0300 

0330 

0400 

0430 

0500 



434 

438 

443 

445 

453 

458 

464 

470 

477 

483 

491 

498 

505 

512 

5 i 9 

526 

532 

539 

545 

552 

558 

564 

570 

576 

581 

587 

592 

597 

601 

606 



278 

279 
279 
278 
276 
274 
271 
268 
267 

267 

268 
271 
273 
276 
278 
281 
283 
285 
287 
289 
291 
293 

295 

296 

298 

299 

300 

301 

302 

303 



8 
8 
9 
9 
9 
9 
9 
10 
10 
11.0 
10 .9 
10 .8 
10.7 
10.6 
10 
10 
10 
10 
10 
10 
10 



4 

4 

3 

0 

1 
0 



9.9 

9.8 

9.5 

9.3 
9 . 0 

8.8 

8.6 
9.0 

9.4 



12 

12 

13 

13 

13 



9 

9 

0 

0 

0 



13.0 
13 . 0 



13 

12 



0 

9 



12 . 9 
13 . 0 
13.0 
13 . 0 
13 . 0 



13 

13 

13 

13 

13 

13 

13 

13 



0 

1 

1 

1 

1 

1 

1 

1 



13. 1 
13. 1 
13. 1 
13 . 1 
13. 1 
13.1 
13. 1 
13.0 



0 

0 

0 



0 

0 



15.0 
15.0 
15.0 
15 
15 . 
15 
15.0 
15.0 
15 . 
15 . 
15.0 
15.0 
15.0 
15.0 
15. 0 
15.0 
15.0 
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Figure 19. Time Evolution of the Sea-Surface Temperature-. 

Case 2. Predicted values are shown by the solid 
line. Observed values are shown by the dashed 
line . 




Figure 20. Time Evolution of Ocean Mixed Laver Depth (MLD) - 
Case 2. * 
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Figure 21. Selected Ocean Vertical Temperature Profiles- 
Case 2. 
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Figure 22. Selected Ocean Vertical Profiles of the 

Horizontal Wind Driven Currents- Case 2. East- 
west component u (solid line) North-south 
ccmponen- v (dashed) . 
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Figure 23. Selected Change in Oceai Vertical Salinicy 
Profiles- Case 2. 
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Figure 26. Time Evolutior. of the Ssa-Surface Temperature-. 

Case 3. Predicted values are show?, by zhe solid 
line. Observed values are shown by tne dashed 
line . 
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Figure 28. Selected Ocean Vertical Temperature Profiles- 
Case 3 . 
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Figure 29. 



Selected Ocean Vertical Profiles of the 
Horizontal Wind Driver. Currents- Case 3. Hast- 
wesT component u (solid line) North-south 
component v (dashed) . 



81 



0 



E 



X 

■»* 

a 



- 10 - 



330 



n7r]inir|rnTpiT[TTr nirp nn] rnTyrm [ m 



5 -20 



-30 



mi 



Prof 1 1 o 1 1 «no« 
900 



1430 

mrpn^ijTirrpTr 






mi 



iinhmliinluiliin 



miliinliiiihit 



-35 0 25 50 75-25 0 25 50 75-25 0 25 50 75 



Salinity Change X 10000 (ppt) 



Fiqure 30. Selected Cha?.ge in Oceaa Vertical Salinity 
Profiles- Case 3. 
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Figure 3 1. Graphical Display of *^he Coaolad Model Predic 
Armospheric Mrxea Layer Quantities- Cass 3. 
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